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O I Hadronic event shape distributions are determined from data in e~^e collisions 

><; ■ between 183 and 207 GeV. From these the strong coupling a^ is extracted in 

^ . 0{al), NLLA and matched (9(a^)+NLLA theory. Hadronisation corrections 

Q-i. evaluated with fragmentation model generators as well as an analytical power 

(— I I ansatz are applied. Comparing these measurements to those obtained at and 
around Mz allows a combined measurement of a^ from all DELPHI data and 

... a test of the energy dependence of the strong couphng. 
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1 Introduction 

Measurements of the strong coupling, a^, of quantum chromodynamics [1] (QCD), the 
theory of strong interaction, using different observables and different analysis methods 
serve as an important consistency test of QCD. Once as is measured at a given scale, 
QCD predicts its energy dependence as described by the renormalisation group equation. 
A measurement of the strong coupling at different scales allows therefore a test of this 
important prediction, which is related to the property of asymptotic freedom [2]. 

At LEP, hadronic final states of the e+e" annihilation are used to study QCD. While 
the process e^e~ -^ qq is described by the electroweak theory alone, the radiation of 
gluons carries sensitivity to properties of the strong interaction. Our analysis uses event 
shape observables to measure the strong coupling. These dimensionless quantities char- 
acterize the topology of the events, e.g. whether the radiation of hard gluons gave rise 
to further jets. 

The strong coupling is measured by comparing experimental cross-sections, ^^, for an 
observable y with the theoretical predictions in which a^ enters as a free parameter. But 
since QCD is the theory of (asymptotically) free quarks and gluons, hadronisation effects 
need to be accounted for. This may be done either with phenomenological models [3-5], 
or with the help of QCD-inspired power corrections [6]. 

The QCD calculation can be performed in different ways as well. The earliest results 
were based on fixed-order perturbation theory [7] . For the observables studied here these 
predictions are limited to the three-jet region. An extension to the four-jet region would 
need next-to-next-to-leading-order (NNLO) corrections to be calculated. On the other 
hand the applicability of these calculations close to the two-jet region is limited as well, 
since in this kinematic domain enhanced logarithms occur [8]. To extend the applicability 
into the two-jet region the summation of these logarithms was developed, the so called 
next-to-leading-log approximation (NLLA) [8,9]. Finally the fixed-order results can be 
combined with NLLA calculations leading to the (9(q;^)+NLLA matched theory for the 
cross-section, R = f ^^dy, [8,9]. According to different "matching schemes", these 
calculations are referred to as e.g. R or logi? matching [8]. 

As a consequence of renormalisation, all perturbative QCD calculations to finite order 
depend upon the renormalisation scale /i, which is an unphysical parameter. The choice of 
H is conventional and the effect of its variation is usually used to estimate the theoretical 
uncertainty. In the NLLA and matched theory even more arbitrary parameters enter, 
related to the phase-space boundary. We will discuss this point and our definition of the 
theoretical uncertainties in section 4.1. 

This paper presents the measurements of event shape distributions in e^e^ collisions 
between 183 and 207 GeV. The data have been reprocessed in 2001 and our final results 
supersede some earlier DELPHI measurements at the corresponding energies [10]. An- 
other change with respect to the previous LEP2 analysis is the use of improved event 
generators for both acceptance correction and background subtraction (see section 2). 

From the event shapes Thrust, C parameter, heavy jet mass, wide and total jet broad- 
ening, as is extracted with four different methods: the differential distributions are com- 
pared to predictions in (9(0;^), pure NLLA and 0{a'^)+NLLA (logR), folded with frag- 
mentation models, and in the fourth method the strong coupling is extracted from the 
mean values using an analytical power correction ansatz. An extension of this analysis 
with respect to the previous one [10] is the use of five observables instead of thrust and 
heavy jet mass only. Additionally the matching procedure for the 0{a1)+NLLA (logR) 
prediction has been modified. For consistency the distributions at Mz from [10] have 



been refitted for these five observables (and with the same fit ranges). The combination 
of five Us values for each method and energy makes the treatment of correlations more 
crucial. Section 4.2 is devoted to this topic. 

From here the analysis proceeds in two steps: first the Us values (together with the 
results from previous measurements at other LEP2 energies and LEPl data) are used to 
test the QCD predicted scale dependence, i.e. to measure the j3 function of the strong 
interaction. Second, assuming the QCD (3 function, all «« values at LEPl and LEP2 
are evolved to a reference energy and combined to a single Ug value for each method. It 
turns out that the weight of the LEP2 data in the combined a^ results is comparable to 
the weight of the LEPl data alone. This unexpected result is due to the fact that at 
LEP2 the bigger statistical uncertainties are compensated for by smaller hadronisation 
and scale uncertainties. 

The paper is organized as follows: in section 2 the selection of hadronic events, the 
determination of the centre-of-mass energy, the correction procedures applied to the data, 
and the suppression of WW and ZZ events are briefly discussed. Section 3 presents event 
shapes and the comparison of the data with predictions from different generators. The 
measurements of Ug from differential distributions are discussed in section 4, while section 
5 describes the as determination from mean values with power corrections. In section 6 
the running of the strong coupling is discussed and section 7 contains the combination of 
all as measurements. Section 8 gives a summary of the results. 

2 Selection and correction of hadronic data 

The analysis is based on data taken with the DELPHI detector in the years from 1997 
to 2000 at centre-of-mass energies between 183 and 207 GeV. Detailed information about 
the design and performance of DELPHI can be found in [11,12]. 

In order to select well-measured charged particle tracks, the cuts given in the upper 
part of Table 1 have been applied. The cuts in the lower part of the table are used to 
select e^e~ -^ Z/7 -^ qq events and to suppress background processes such as two- 
photon interactions, beam-gas and beam-wall interactions, leptonic final states, events 
with hard initial-state radiation (ISR), WW and ZZ pair production. 

At energies well above Mz the high cross-section of the Z resonance raises the proba- 
bility of events with hard ISR. These "radiative return events" constitute a large fraction 
of all hadronic events. The initial-state photons are typically aligned along the beam 
direction and are identified inside the detector only at a rate of about 10% . In order 
to evaluate the effective hadronic centre-of-mass energy of an event, considering ISR, an 
algorithm called Sprime is used [13]. Sprime is based on a 3C fit imposing transverse 
momentum and energy conservation. Several assumptions about the event topology are 
tested. The decision is taken according to the x^ obtained from the constrained fits with 
different topologies. 

Figure l(left) shows the spectra of the calculated energies for simulated and measured 
events after all but the ys' cut. A cut on the reconstructed centre-of-mass energy ys' > 
90%y^ is applied to discard radiative return events. 

Two-photon events are strongly suppressed by the cuts. Leptonic background was 
found to be negligible in this analysis as well. 

Since the topological signatures of QCD four-jet events and hadronic WW, ZZ and 
other events with four-fermions (4F) in the final state are similar, no highly efficient 
separation of QCD events and backgrounds is possible. Furthermore any 4F rejection 



Track 


0.2GeV/c<p<100GeV/c 


selection 


Ap/p < 1.0 




measured track length > 30 cm 




distance to LP in r0 plane < 4 cm 




distance to LP. in 2; < 10 cm 


Event 


-'''charged ^ ' 


selection 


25° < ^Thrust < 155° 




Etot > 0.50v/i 




V? > 90%v^ 




A^chargcd>5005^in+1.5 




-' 'charged ^ 4z 



Table 1: Selection of tracks and events, p is the momentum, Ap its error, r the radial 
distance to the beam-axis, z the distance to the beam interaction point (LP.) along the 
beam-axis, the azimuthal angle, A^charged the number of charged particles, 6'Thrust the 
polar angle of the thrust axis with respect to the beam, i?tot the total energy carried 
by charged and neutral particles, vi' the reconstructed centre-of-mass energy, ^/s the 
nominal centre-of-mass energy, and -Bmin is the minimal jet broadening. The first two 
cuts apply to charged and neutral particles, while the other track selection cuts apply 
only to charged particles. 



implies a severe bias to the QCD event shape distributions, which needs to be corrected 
by simulation. 

Our suppression of these backgrounds uses a two-dimensional cut in the plane 
spanned by the charged particle multiplicity (A^chargcd) and the narrow jet broadening 
-Bmin = min(i?_|_, i?_). B± is defined as the normalized sum over the transverse momen- 
tum of charged and neutral particles in the two event hemispheres separated by the plane 
perpendicular to the thrust axis tit'- 



B. 



^ \pi X Ut 




2Eip"^ 



By applying a cut on an observable calculated from the narrow event hemisphere only, the 
bias to event shape observables mainly sensitive to the wide event hemisphere is reduced. 
The charged particle multiplicity is used to reduce the 4F contribution further. The two- 
dimensional cut in the A^charged--Bmin plauc cxploits the different correlation between these 
observables for QCD and four-fermion events, as shown in Figure 1 (right). Especially 
some reduction for semi-leptonic decaying 4F events is gained. The lines indicate the 
cut values chosen. This cut suppress almost 90% of the four-fermion background. The 
remaining 4F contribution is estimated by the WPHACT [14] generator and subtracted 
from the measurement. 

Table 2 contains the integrated luminosities at different energies, the cross-sections 
for signal and background and summarizes the selection statistics. The cross-sections 
were taken from the simulation which was used to correct the data and to subtract the 
background. The cross-sections for the 4F background are quoted for charged current 



V^ [ GoV] 
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189 


192 
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205 


207 


i:[pb]-' 


55.73 


157.97 


25.34 


67.29 


78.07 


39.31 


76.33 


130.12 


-JiV"[pb] 


108.78 


100.05 


96.06 


91.31 
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79.78 


-T.^o^Sp^] 


23.09 


21.24 
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19.36 
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a^^'--'-'[pb] 


17.54 
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19.10 
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19.97 


20.10 


20.14 


a^^'^^'lph] 


8.16 


8.15 


8.14 


8.08 


8.03 


8.01 


7.93 


7.90 


<^HE 


0.721 


0.720 


0.736 


0.740 
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0.734 


0.736 


0.749 


<^CC 


0.090 


0.100 
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0.124 


<^NC 


0.017 


0.017 


0.017 


0.017 


0.015 


0.016 


0.016 


0.016 


Phe,qcd 


0.867 


0.848 


0.837 


0.828 


0.808 


0.801 


0.790 


0.795 


^HE ■ PHE,QCD 


0.625 


0.610 


0.617 


0.612 


0.594 


0.588 


0.581 


0.593 


# selected events 


1070 


2848 


455 


1164 


1303 


653 


1203 


2036 


# CC background 


87.8 


296.0 


50.1 


147.7 


189.8 


94.1 


195.6 


315.4 


# NC background 


7.7 


21.4 


3.54 


9.36 


9.65 


4.88 


9.67 


16.60 



Table 2: Luminosities, cross-sections of QCD signal and background from four-fermion 
events (split into neutral current, NC, and charged current, CC), selection efficiencies, 
e, and purities, p. The subscript HE denotes QCD high energy events, i.e. with 
vs' > 0.9 ■ ^/s. Also given is the total number of selected events and the expected 
number of remaining four-fermion events. 



(CC) and neutral current (NC) contributions separately. Details on the four-fermion 
simulation in DELPHI can be found in [15]. 

The influence of detector effects was studied by passing events (generated with /C/C 
[16]) and fragmented with Jetset/Pythia [3] using the DELPHI tuning described in [17] 
through a full detector simulation (Delsim [11]). This simulation is improved with 
respect to the previous LEP2 analysis [10] by including electroweak corrections (multiple 
photon emission, treatment of ISR and FSR etc.). These simulated events are processed 
with the same reconstruction program and selection cuts as are the real data. In order to 
correct for cuts, detector, and ISR effects a bin-by-bin acceptance correction C, obtained 
from e~^e~ -^ Z/7 -^ qq simulation, is applied to the data: 



a 



h{f. 



J/gcn,noISR 



Kfi) 



(2) 



where /i(/i)gen,noiSR represents bin i of the shape distribution / generated with the tuned 
generator. The subscript noISR indicates that only events without relevant ISR {^/s — 
V^ < 0.1 GeV) enter the distribution. /i(/j)acc represents the accepted distribution / as 
obtained with the full detector simulation. The more detailed matrix correction used for 
the data measured at the Z peak [18] is not applied here, because of the smaller statistics 
at LEP2. 



3 Event shape distributions and mean values 

Selected event shape distributions at 189 and 207 GeV are shown in Figures 2 and 3. 
The definitions of these observables are given is Section 4. 

The data in Figures 2 and 3 are corrected to be comparable with e~^e~ -^ Z/7 -^ qq 
simulation of charged and neutral hadron production. In the data all charged particles are 
assumed to have pion mass while neutral particles are considered massless. The monte 
carlo correction of the data includes the effects of the simulated particle masses. The 



DELPHI 196 GeV 

• data 
QCD+ WW/ZZ MC 



rt 



^^^ 





- N^ = SOOB .+1.5 

- N -7 

eV -> WW -* hadrons + Iv 



u c" -* WW — > luiilroiis 



eV — > Z/y — * hadrons 



lOU 120 140 160 180 200 220 



0.02 0.04 dM 0.0S 0.1 0.12 0.14 0.16 0.18 0.2 



Vs' [GeV] 



Figure 1: Left: reconstructed centre-of-mass energy vs' . Right: simulation of 
four-fermion background and QCD events in the A'"chargcd--Bmin plane. The lines delin- 
eate the accepted region. 



Figures 2 and 3 compare these data with the Jetset [3], Ariadne [4] and Herwig [5] 
generators as tuned by DELPHI [17] with LEPl data. The amount of 4F-background 
which was subtracted to obtain the final data points is also shown. The acceptance 
corrections are plotted in the upper inset. 

The Tables 8 and 9 at the end of the paper contain mean values and higher moments 
for the event shapes 1-T, C parameter, Ml/E^^^, -Bmax and Bsum- Also included are the 
results for alternative definitions of the heavy jet mass as proposed in [19]. They are 
obtained if in the definition of the heavy jet mass the invariant mass is calculated with 
the following replacements: 



or {Ei,pi 



{\Pi\,Pi) 



In what follows we will refer to these observables as p-scheme and E-scheme definitions 
of the heavy jet mass. 

In order to estimate the systematic uncertainty from the selection and correction 
procedure, the effects of the following changes with respect to the standard values have 
been considered: Ach±l, 0thrust=t5° and \^/^/s±0.025. For the 4F cross-section a change 
of ±5% has been considered and the uncertainty due to the acceptance correction was 
estimated by a change of ±0.02. The last uncertainty agrees with that of [17] where the 
systematic uncertainty was verified using independent data. Half of the difference between 
up- and downward variation is regarded as one component of the systematic uncertainty. 
These five contributions are added in quadrature to estimate the experimental systematic 
uncertainty. 





/ \ 




,, _^^ 


y 


(' ''~~~- 


— ' ~-^- 






^ 









-'d) 




DELPHI 
189 GeV 




r ^r^ 


"^t 


• data 
— Jetset7.4PS 




'* 


-v 


HerwigS.S 








^>s^^^ Ariadne 4.08 




■ fi 

" '1 

■ vi 

- :i 

■ ■! 




^^^ \ 


x* 


■ ■'! 




r 




\'> 











Figure 2: Event shape distributions of 1-Thrust (1 — T), heavy jet mass [M'^jE'l-^^^ wide 
jet broadening (-Bmax) and total jet broadening (i?sum) at 189 GeV. The upper inset shows 
the acceptance corrections. The central part shows data with statistical uncertainties, 
simulation and the four-fermion background which was subtracted from the data. 
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Figure 3: Event shape distributions of 1-Thrust (1 — T), heavy jet mass (M^/i^^jg), wide 
jet broadening (i?max) and total jet broadening (i?sum) at 207 GeV. The upper inset shows 
the acceptance corrections. The central part shows data with statistical uncertainties, 
simulation and the four-fermion background which was subtracted from the data. 
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Figure 4: Fit ranges for the different observables and methods to determine ag- 

4 Determination of ols from event shape distribu- 
tions 

Our determination of as is based on the five variables 1-Thrust (1 — T), C-parameter, 
heavy jet mass (M^/_E^jg), wide jet broadening (-Bmax) and total jet broadening (-Bsum)- 
The thrust, T, is defined as: 

T=max(^-'^-:"'U^if-:^^' . 

" I Y.i\Pi\ J Y.i\Pi\ 

The vector which maximizes the above expression defines the thrust axis, n-y. The plane 
perpendicular to the thrust axis divides the event into two hemispheres. Based on this 
separation several other event shapes can be defined. One defines the heavy jet mass by 
the following expression: 



Mh7E4 = max(M7M^)/S, 
M^ denotes the invariant mass of the two hemispheres: 

^±Pi-nT>0 

Here Pi is the four-momentum of the ith particle. Using the expression B± as defined in 
Equ.l the wide and total jet broadening are defined as: 

5max = max(5+,B_) 

-Bsum = -B+ + -B- 

The linear momentum tensor offers a possibility to define event shapes without dis- 
tinguishing an event axis. It is defined as: 

e'^'=^^/^bl| with: a,b = x,y,z 



From its eigenvalues Aj the C-parameter is defined: 

C = 3(AiA2 + A1A3 + A2A3) 

From these differential distributions as is determined by fitting an a^-dependent QCD 
prediction folded with a hadronisation correction to the data. The following QCD pre- 
dictions are used: (9(a^), pure NLLA, and the modified (9(q;^)+NLLA in the logi?- 
scheme [7-9,20,21]. Hadronisation corrections are calculated using the Jetset PS model 
(Version 7.4 as tuned by DELPHI [17]). In each bin the QCD prediction is multiplied by 
the hadronisation correction 



a 



had 



f Sim. 
ihad 
^Sim. 
./part 



(3) 



where /^a™' (/part') is the model prediction on hadron (parton) level. The parton level is 
defined as the final state of the parton shower created by the simulation. 

The fit ranges used for the different QCD predictions are shown in Figure 4. The 
lower edges are chosen in such a way, that the hadronisation corrections in the 2-jet 
region remain small (< 10%) for LEP2 energies. The upper limit of the fit ranges ensures 
that the signal-to-background ratio is above 1. The ranges for pure NLLA and 0{a1) fits 
are chosen to be distinct, so that the results are statistically uncorrelated. 

In [18] it has been shown that fixing the renormalisation scale to /x = ^/s results in a 
poor description of the data. Therefore, the experimentally optimized scales (fiEos) 



observable 


experimentally optimized 




scales {Xfj_ = fi/y/s) 


1-T 


0.057 


C 


0.082 


MhVi?v\ 


0.060 


-Dmax 


0.143 


-Dsum 


0.096 



Table 3: The values for the experimentally optimized scales from [18]. 

from [18] (see Table 3) are used for the 0{a'j.) fits. For the NLLA and the combined 
NLLA+0{a'^) fits, n is still set equal to y/s. This is the conventional choice of scale 
for resummed and matched calculations and allows a direct comparison with the results 
from other experiments [22]. Furthermore the meaning of the renormalisation scale n in 
resummed calculations is different from its interpretation in the framework of fixed-order 
perturbation theory. While in Ola^) fi paramatrizes the choice of the renormalisation 
scheme this interpretation is lost in resummed calculations. This difference makes the 
application of experimentally optimized scales especially for NLLA+0{a'^) predictions 
meaningless. Tables 10-18 at the end of this paper contains all as values derived from 
event shapes. 



4.1 Definition of uncertainties 

Experimental systematic uncertainties are obtained from fits to distributions evalu- 
ated with different cuts and corrections. These variations are described in section 3. 
The hadronisation uncertainty is taken to be the bigger of the two differences when the 
hadronisation correction is determined from the Ariadne [4] and Herwig models [5] 
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alternatively. The Jetset result is used as the central value. In all cases the dominant 
systematics come from the theoretical uncertainty. The conventional method for estimat- 
ing this uncertainty is to consider the effect of a renormalisation scale variation [10] when 
fitting the experimental distributions. This method, however, has at least two draw- 
backs: since the resulting scale uncertainty is positively correlated with the measured 
as, this definition produces a bias towards small a^ values when combining the results 
of e.g. different observables. Secondly there are indications that observables calculated 
only in one hemisphere (like the heavy jet mass or -Bmax) yield less reliable results in 
the resummation of leading logarithms [23]. This should be reflected in their theoretical 
uncertainty. Conversely the scale variation yields the smallest uncertainty for the heavy 
jet mass and especially -Bmax- For these reasons a new definition of the theoretical uncer- 
tainty for the logR prediction was developed in cooperation with the LEP QCD working 
group. By construction, the NLLA calculations do not vanish at the phase-space limit 
Umax [8]. In the so-called modified theory (NLLA or matched) they are forced to vanish 
by the replacement: 

L = In - ^ L = In 

y 



X -y X -y^ 

In agreement with the LEP QCD working group ymax is chosen as the maximum value 
of the parton shower simulation [24]. Usually X = 1 is chosen for the quantity X, as 
suggested by the authors of [8], although different values for this X scale introduce only 
subleading contributions [25] . The theoretical uncertainty of the logR prediction in this 
analysis is now defined as half of the difference when X is varied between 2/3 and 3/2. 
By this new definition of the uncertainty the observables M^/E'^-^^ and -Bmax, which are 
calculated in one hemisphere only, get a bigger uncertainty compared to the uncertainty 
estimated by fi variation. The same definition of the theoretical uncertainty has been 
adopted for the pure NLLA prediction. 

For the C(q;s) calculation we use, as in the previous publication [10], the effect from 
the variation around the experimentally optimized scales, ^eoS) between Q.b^Eos and 
2iiEos to estimate the theoretical uncertainty. 

In order to avoid the effect mentioned above of a positive correlation, all scale varia- 
tions have been calculated for a fixed value of as from the theoretical distributions for 
each method separately. The fixed as value is chosen as the average as value of the 
combination. To obtain this value the procedure has to be iterated. 

4.2 Method for combining the cxg measurements 

For a combination of the as results from different observables calculated from the same 
data sets a proper treatment of the correlation is mandatory. The average value y for 
correlated measurements y, is [26]: 

y = 2_^ Wiyi with: Wi 



Y.uAv~-'\ 



kl 



Note that the weights w can be negative, if the correlation pij between two quantities 
i and j is bigger than Oijoj. Here a is the uncertainty of the corresponding quantity 
with Gi > Gj. The covariance matrix V has an additive structure for each source of 
uncertainty: 

T7 T^stat I T^sys.cxp. _i_ T/had _i_ T^scalc 
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Its statistical component is estimated with simulation which yields correlations of typ- 
ically > 80%. The correlation of systematic uncertainties is modeled by the minimum 
overlap assumption: 

Vij = mm{ai,af). 

The as values evaluated from the distributions and their mean values taking correlations 
into account are given in the Tables 10-18 at the end of the paper. 

5 Determination of cXs from mean values with power 
corrections 

The analytical power ansatz for non-perturbative corrections by Dokshitzer and Web- 
ber [6,27] including the Milan factor established by Dokshitzer et al. [28,29] is used to 
determine as from mean event shapes. This ansatz provides an additive term to the 
perturbative 0{af) QCD prediction: 

(/) = —/ f^da = if,,,,) + (/po.) , (4) 

o"tot J do- 

where the 2nd order perturbative prediction can be written as 



2n \ s ' 7 V 27r 

with A and B being the perturbative coefficients [7,30] , n being the renormalisation scale 
and bo = (33 — 2Nf)/12TT. The power correction is given by 



(/pow) = Cf—^M-^^ 



ao{iJ,i) - asifi) - 6o ■ log^ + 7— + 26o a^(/i) 






where ao is a non-perturbative parameter accounting for the contributions to the event 
shape below an infrared matching scale fij and K = (67/18 — tt^/6)Ca — 5Nf/9. The 
Milan factor A^ is set to 1.49, which corresponds to three active flavours in the non- 
perturbative region. The observable-dependent quantities A, B and Cf are listed in Table 
4. For the jet broadenings C/ takes a more complicated form [31]: 

Here cb is 0.5 or 1 for (-Bmax) or (-Bsum) respectively, 770 = —0.6137. The infrared 

matching scale is set to 2 GeV as suggested by the authors of [6] , the renormalisation 
scale /i is set to ^ i.e. the MS scheme is used, since the power corrections are provided 
only in this scheme. 

Besides as these formulae contain ao as the only free parameter. In order to measure as 
from the high energy data this quantity has to be determined. To infer ao, a combined 
ffi of as and ao to a large set of measurements at different energies [32] is performed. 
For ^/s > Mz only DELPHI measurements are included in the ffi. Figure 5 (left) shows 
the measured mean values of our five observables as a function of the centre-of-mass 
energy together with the results of the fit. The resulting values of ao arc summarized in 
Table 5. The ffist uncertainty in Table 5 is taken from the ffi to the data with full errors. 
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observable 


^/ 


Bf 


c/ 


(1-T) 


2.103 


44.99 


2 


(C) 


8.638 


146.8 


37r 


{Ml/El,) 


2.103 


23.24 


1 


\-Dmax/ 


4.066 


-9.53 


Eq.5 


(-Bsum) 


4.066 


64.24 


Eq.5 



Table 4: A and B coefficients for the expansion of the mean values in Q;s/27r, and values 
for the observable dependent c/. 

while the second uncertainty reflects the effect of a variation 0.5/i < /i < 2/i. Figure 5 
(right) shows the flt results also in the cts-ao plane. The extracted a^ values are supposed 
to be observable independent and around 0.5 [27,29]. However, higher order effects are 
expected to violate this universality. Within the theoretically expected accuracy of 20% 
this universality is fulfllled. 

After flxing ao for each observable to the values in Table 5, the as values correspond- 
ing to the high energy data points can be calculated from Eq. (4). The effect of an a^ 
variation within its uncertainty was found to be well within the systematic uncertainties 
of as- By using the a^ value from the global flt, the determination of a^ uses the DELPHI 
data points twice. But since the global flt is dominated by the low-energy data the effect 
is negligible, as is calculated for all observables individually and then combined taking 
correlations into account as described in section 4.2. An additional scale uncertainty is 
calculated by varying /i for a flxed value of as and the infrared matching scale /x/ from 
1 GeV to 3 GeV. The as results are summarized in the Tables 19 to 22 at the end of 
the paper. The total error for this method is smaller than e. g. for NLLA+(9(q;^) flts. 
However, the hadron level which is experimentally accessible does include the effects of 
resonance decays and hadron masses which are not accounted for in the calculation of 
power corrections. In order to investigate the influence of different hadron level deflni- 
tions a Monte Carlo study was performed in [19]. Three different hadron level deflnitions 
were considered: (i) hadrons which are primary produced, (ii) stable hadrons after reso- 
nance decays or (iii) particles out of a subsequent decay into two massless particles. The 
subsequent determination of the strong coupling from power corrections leads to a shift 
in as of ±0.0035 with respect to the hadron level (ii). With regard to these extreme 



Observable 


ao(2GcV) 


as{Mz) 


xVndf 


(1-T) 


0.532 ± 0.011 ± 0.002 


0.122 ± 0.001 ± 0.009 


69/43 


{C) 


0.442 ± 0.010 ± 0.008 


0.126 ± 0.002 ± 0.006 


18/22 


{Ml/El,) 


0.620 ± 0.028 ± 0.010 


0.119 ± 0.002 ± 0.004 


10/32 


{Ml/ El,) (E def) 


0.576 ± 0.113 ± 0.002 


0.111 ± 0.005 ± 0.003 


5/14 


{Ml/ El,) (p def) 


0.517 ±0.110 ±0.003 


0.110 ± 0.005 ± 0.004 


3/14 


X-t'max/ 


0.460 ± 0.029 ± 0.078 


0.116 ± 0.001 ± 0.002 


7/22 


(-Bsum) 


0.452 ± 0.014 ± 0.015 


0.118 ± 0.001 ± 0.004 


12/22 



Table 5: ao and as values from the global flt of the Dokshitzer- Webber ansatz for mean 
values to e+e^ data from several experiments [32]. Only the ao values are used further 
for the as determination from single mean values at LEP2. 
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Figure 5: Left: Dokshitzer- Webber fit to several mean values. The dotted line shows the 
perturbative contribution. Right: the results of the global fits in the as-ao plane. The 
vertical line with shading shows the world average of as- 



assumptions one may therefore assign 0.007/v 12 = 0.002 as an additional uncertainty 
which accounts for the fact that resonance decays and hadron masses are not considered 
in the calculation. 



6 The running of cXs 

The as values determined at different energies are used to test the predicted scale 
dependence of the coupling. We include also the LEP2 results at 133, 161 and 172 GeV 
from [10]. For a^ at and around Mz we have reanalyzed the distributions from [10] for 
the five observables and combined the results using the same treatment of correlations 
as described in section 4.2. For as from mean values the measurements of events with 
reduced centre-of-mass energy between 44 and 76 GeV [33] and the data between 133 and 
172 GeV [10] have been included as well. In the Tables 10-22 at the end of this paper all 
these as values are provided. 

The logarithmic energy slope of the inverse coupling is given by: 



da„ 



with bo 



33-2N 



and bi 



dlogi/s 

153-19Af/ 



2bo + 2bias + 



(6) 



corresponding to the first coefficients of the (3 function. 



127r """"" ""-L 24?;^ 

The measurement of this quantity allows both a test of QCD and a consistency check 
of the four different methods used to determine as- Equation 6 shows that in leading 
order da^^/dlog y/s is independent of a^ and twice the first coefficient of the /3 function. 
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Evaluating this equation in second order results in a slight dependence on a^. With 
Q;s=0.11 (which corresponds to Aqcd = 230 MeV and y/s =150 GeV, the average energy 
of our measurements) one obtains da~^/d\og^/s = 1.27. 

Table 6 gives the slopes when fitting the function {blog^/s + c) to the a~^ values. 
The correlation between the as measurements is taken into account by including the full 
covariance matrix in the definition of the x^ function. The correlation is modeled as 
described in section 4.2. The only difference here is that the statistical uncertainties are 
uncorrelated. The a^ values and the fit of their energy dependence are also displayed 
in Figure 6. The results are in good agreement with the QCD expectation. Using the 
definition of the bi the result for the slope can be converted into the number of active 
flavours, Nf. These numbers are also included in Table 6. 

A model-independent way to measure the [3 function is offered by applying the renor- 
malisation group invariant (RGI) perturbation theory to the mean values of event shapes 
directly 
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Figure 6: Energy dependence of as as obtained from event shape distributions using 
different theoretical calculations. The total and statistical (inner error-bars) uncertainties 
are shown. The band displays the average values of these measurements when extrapo- 
lated according to the QCD prediction. The dashed lines show the result of the 1/ log y/s 
fit. 
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Figure 7: Results of combining all DELPHI Ug measurements at LEPl and LEP2. 
The total and statistical (inner error-bars) uncertainties of the individual measurements 
are displayed. The central results are the correlated means. For comparison also the 
unweighted and total-error weighted (but uncorrelated) averages are shown. For the un- 
weighted mean the size of the error-bars indicate the RMS of the measurements. The 
weights of the individual measurements within the correlated average are given in brack- 
ets. Note that these can turn negative in the presence of strong correlation as e.g. for 
the a» results from mean values. 
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theory used for measurement 


dt^v-s ± ^t^t ± 'y' 


x'/ndf 


Nf 


mean values + power corr. 
C(a2)+NLLA (logR) 

NLLA 


1.11 ± 0.09 ± 0.19 
1.32 ± 0.11 ± 0.27 
1.27 ±0.15 ±0.33 
1.40 ± 0.17 ± 0.44 


1.25 
0.58 
0.29 
0.83 


6.3±1.7 
4.6±2.3 
5.0±2.9 
4.0±3.8 


QCD expectation 


1.27 




5 



Table 6: Results for the slope b when fitting the function l/{b\ogy/s + c) to as values 
obtained for the different energies. Also given is the corresponding result for the number 
of active flavours, Np- 

7 Combination of all DELPHI cXg measurements 

As shown in the last section, the energy dependence of as is shown to be in good 
agreement with the QCD prediction. Assuming now the validity of QCD, all as results 
can be evolved to a reference energy, e.g. Mz, and combined to a single as{Mz) mea- 
surement. Again we include results from other LEP2 energies and LEPl as described in 
section 6. Combining the as results is, again, complicated by correlations among these 
measurements. Although the measurements at different energies are clearly statistically 
independent, the systematic and theoretical uncertainties are not. Again this part of the 
covariance matrix was modeled assuming minimum overlap. 

The results of the combinations are given in Table 7 and displayed in Figure 7. The 
Figure contains in brackets also the weights of the individual measurements within the 
average. As can be seen from these numbers the weight of LEPl and LEP2 measurements 
are roughly the same, since smaller theoretical uncertainties at LEP2 compensate for the 
larger statistical error. As can be seen from Table 7 the total error is still dominated by 
the scale uncertainty. The result with the smallest total uncertainty is deduced from the 
0{a1) prediction from distributions. Here the total uncertainty is 0.0033. This value 
can be compared with the central result of the DELPHI analysis [18] from the observable 
jet cone energy fraction (JCEF) alone: q;s=0.1180±0.0018. The different precision is 
mainly due to the definition of the scale uncertainty. While we use the variation 0.5/i < 
jj, < 2yU, the analysis [18] changes the corresponding quantity only between V0.5/i and 
v2/i. The actual choice of the /i variation is not guided by solid theoretical arguments 



theory 


as{Mz) 


stat. 


sys.exp. 


had. 


scale 


tot 


mean values ± power corr. 


0.1184 


0.0004 


0.0008 


0.0022 


0.0031 


0.0039 


C(a2)+NLLA (logR) 


0.1205 


0.0010 


0.0018 


0.0013 


0.0048 


0.0054 


0{al) 


0.1157 


0.0008 


0.0016 


0.0016 


0.0022 


0.0033 


NLLA 


0.1093 


0.0012 


0.0020 


0.0011 


0.0050 


0.0056 



Table 7: Results of combining all DELPHI as measurements at LEPl and LEP2. For the 
as results from mean values with power corrections "hadronisation uncertainty" (had.) 
denotes the combined effect of the /i/ variation and the uncertainty related to resonance 
decays and particle masses as described at the end of Section 5. The scale uncertainty is 
either the effect of a variation of the renormalisation scale /x {0{a'^) and power corrections) 
or the effect of changing the X scale (see Section 4.1). 
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and was studied in all possible detail in [18]. However, our definition of the theoretical 
uncertainty yields good agreement with the average root-mean-square of the fits to the 
five different observables at the same energy (see Tables 10-12). Another reason for 
the higher precision in [18] is the use of the observable JCEF, which has particularly 
small uncertainties from hadronisation and scale variation. However, the focus of this 
work is the analysis of five observables with several different techniques. At present the 
theoretical uncertainties of a^ measurements from event shape distributions are subject 
of a debate. Further substantial progress can only be achieved by the arrival of next-to- 
next-to-leading-order (NNLO) calculations. 

8 Conclusion 

A measurement of event shape distributions and mean values is presented as obtained 
from data at centre-of-mass energies from 183 to 207 GeV. The strong coupling constant 
Us has been determined from the event shape variables Thrust, C parameter, heavy 
jet mass, wide and total jet broadening, with four different methods: the differential 
distributions are compared to predictions in 0{a1), pure NLLA and 0(q;^)-|-NLLA (logR), 
folded with fragmentation models, while from the mean values, as is extracted using an 
analytical power correction ansatz. The Ug values are combined with results obtained at 
other LEP2 energies and at and around Mz- This allows both a combined measurement 
of Us and a test of the running of Ug- In these combinations the full correlation between 
energies and observables was taken into account. 

The main aim of this study is the comparison of different methods to extract q;^ and 
its scale dependence. Within their uncertainties all techniques yield consistent results. 
The Us with smallest uncertainty is obtained from 0{a1) with experimentally optimised 
scales: 

as{Mz) = 0.1157 ± 0.0008 (stat) ± 0.0016 (sys.ex.) ± 0.0016 (had) ± 0.0022 (scale) 
= 0.1157 ±0.0033 (tot) 

The current world average from the particle data group is 0.1172 ± 0.0020 [34]. 

For the energy dependence of the strong coupling the highest precision is obtained for 
the as values derived from mean values with power corrections: 

"" 1.11 ±0.09 (stat) ± 0.19 (sys) 



dlogv^s 
The last number has to be compared with the QCD expectation of 1.27. 
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V^ 


(l-T) 


{{1-Tf) 


{(i-r)«) 


183 
189 
192 
196 
200 
202 
205 
207 


0.0592± 0.0024± 0.0020 
0.0557± 0.0016± 0.0022 
0.0502± 0.0040± 0.0023 
0.0592± 0.0029± 0.0024 
0.0541± 0.0028± 0.0025 
0.0480± 0.0040± 0.0025 
0.0446± 0.0030± 0.0026 
0.0536± 0.0023± 0.0027 


0.00766±0.00070±0.00054 
0.00658±0.00048±0.00060 
0.00454±0.00116±0.00064 
0.00810±0.00085±0.00067 
0.00613±0.00086±0.00071 
0.00330±0.00121±0.00072 
0.00322±0.00090±0.00077 
0.00572±0.00066±0.00079 


0.00154 ±0.00021±0.00015 
0.00121 ±0.00015±0.00017 
0.00055 ±0.00035±0.00018 
0.00171 ±0.00026±0.00019 
0.00101 ±0.00027±0.00020 
0.00003 ±0.00039±0.00020 
0.00019 ±0.00027±0.00022 
0.00086 ±0.00020±0.00022 




(C) 


iicr) 


{(cr) 


183 
189 
192 
196 
200 
202 
205 
207 


0.2286± 0.0070± 0.0106 
0.2304± 0.0046± 0.0113 
0.2060± 0.0115± 0.0117 
0.2181± 0.0080± 0.0121 
0.2139± 0.0079± 0.0126 
0.2066± O.Ollli 0.0127 
0.1726± 0.0088± 0.0133 
0.2081± 0.0065± 0.0136 


0.08846±0.00544±0.00952 
0.09252±0.00369±0.01030 
0.06634±0.00915±0.01073 
0.08091±0.00657±0.01118 
0.07882±0.00658±0.01170 
0.06730±0.00919±0.01184 
0.03792±0.00743±0.01251 
0.07123±0.00540±0.01278 


0.04585 ±0.00413±0.00824 
0.05114 ±0.00287±0.00894 
0.02676 ±0.00716±0.00933 
0.03909 ±0.00515±0.00973 
0.03907 ±0.00524±0.01020 
0.02761 ±0.00737±0.01032 
0.00349 ±0.00596±0.01092 
0.03150 ±0.00426±0.01116 




(-Bsum) 


{(B=unO'> 


{(B^umf) 


183 
189 
192 
196 
200 
202 
205 
207 


0.0953± 0.0023± 0.0010 
0.0920± 0.0015± 0.0010 
0.0893± 0.0038± 0.0010 
0.0931± 0.0026± 0.0010 
0.0927± 0.0026± 0.0010 
0.0954± 0.0035± 0.0010 
0.0845± 0.0028± 0.0010 
0.0902± 0.0021± 0.0010 


0.01334±0.00070±0.00028 
0.01192±0.00047±0.00030 
0.01113±0.00117±0.00031 
0.01266±0.00082±0.00032 
0.01254±0.00081±0.00034 
0.01344±0.00111±0.00034 
0.00952±0.00086±0.00036 
0.01151±0.00066±0.00036 


0.00247 ±0.00020±0.00011 
0.00199 ±0.00013±0.00012 
0.00178 ±0.00034±0.00013 
0.00224 ±0.00024±0.00014 
0.00222 ±0.00023±0.00015 
0.00257 ±0.00032±0.00015 
0.00131 ±0.00025±0.00016 
0.00188 ±0.00019±0.00016 




(Bmax) 


((Bo.a.)^) 


((Bmax)^) 


183 
189 
192 
196 
200 
202 
205 
207 


0.0663± 0.0021± 0.0021 
0.0652± 0.0014± 0.0022 
0.0621± 0.0035± 0.0022 
0.0668± 0.0024± 0.0022 
0.0659± 0.0024± 0.0023 
0.0666± 0.0033± 0.0023 
0.0625± 0.0025± 0.0023 
0.0658± 0.0020± 0.0023 


0.00688±0.00053±0.00034 
0.00656±0.00037±0.00036 
0.00557±0.00096±0.00038 
0.00719±0.00064±0.00039 
0.00699±0.00063±0.00041 
0.00671±0.00089±0.00041 
0.00585±0.00068±0.00044 
0.00695±0.00053±0.00044 


0.00095 ±0.00013±0.00007 
0.00089 ±0.00009±0.00008 
0.00061 ±0.00026±0.00008 
0.00105 ±0.00016±0.00009 
0.00100 ±0.00016±0.00009 
0.00087 ±0.00023±0.00009 
0.00073 ±0.00017±0.00010 
0.00100 ±0.00013±0.00010 



Table 8: Mean values and higher moments of the Thrust, C, B 
with statistical and systematic errors. 



max and Bsum distributions 
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V^ 


(MhVSv'i.) 


{(MiVS^iJ^) 


{(M,VS?i.)-') 


183 
189 
192 
196 
200 
202 
205 
207 


0.0457± 0.0023± 0.0012 
0.0437± 0.0016± 0.0013 
0.0406± 0.0039± 0.0013 
0.0441± 0.0027± 0.0014 
0.0451± 0.0027± 0.0015 
0.0460± 0.0038± 0.0015 
0.0401± 0.0028± 0.0016 
0.0444± 0.0022± 0.0016 


0.00451±0.00066±0.00027 
0.00408±0.00045±0.00030 
0.00285±0.00117±0.00032 
0.00421±0.00079±0.00034 
0.00458±0.00078±0.00036 
0.00470±0.00112±0.00036 
0.00338±0.00080±0.00039 
0.00439±0.00066±0.00040 


0.00068 ±0.00021±0.00006 
0.00060 ±0.00014±0.00007 
0.00024 ±0.00040±0.00008 
0.00060 ±0.00024±0.00008 
0.00071 ±0.00025±0.00009 
0.00083 ±0.00037±0.00009 
0.00045 ±0.00023±0.00009 
0.00068 ±0.00021±0.00010 




{Mi/E^^,,p) 


{{Mi/E'i,,pf) 


((M,?/S2j^p)«) 


183 
189 
192 
196 
200 
202 
205 
207 


0.0427± 0.0023± 0.0012 
0.0411± 0.0016± 0.0013 
0.0384± 0.0039± 0.0013 
0.0413± 0.0027± 0.0014 
0.0424± 0.0027± 0.0015 
0.0436± 0.0038± 0.0015 
0.0380± 0.0029± 0.0016 
0.0420± 0.0023± 0.0016 


0.00421±0.00066±0.00027 
0.00383±0.00045±0.00030 
0.00274±0.00117±0.00032 
0.00396±0.00079±0.00034 
0.00426±0.00078±0.00036 
0.00451±0.00112±0.00036 
0.00320±0.00080±0.00039 
0.00419±0.00067±0.00040 


0.00068 ±0.00021±0.00006 
0.00060 ±0.00014±0.00007 
0.00025 ±0.00039±0.00008 
0.00057 ±0.00024±0.00008 
0.00065 ±0.00025±0.00009 
0.00081 ±0.00037±0.00009 
0.00043 ±0.00023±0.00009 
0.00065 ±0.00021±0.00010 




{Ml/El,^E) 


{{M?jE^^,,Er) 


{{M^/E'i,,Er) 


183 
189 
192 
196 
200 
202 
205 
207 


0.0434± 0.0023± 0.0012 
0.0417± 0.0016± 0.0013 
0.0391± 0.0039± 0.0013 
0.0420± 0.0027± 0.0014 
0.0430± 0.0027± 0.0015 
0.0440± 0.0038± 0.0015 
0.0385± 0.0029± 0.0016 
0.0426± 0.0023± 0.0016 


0.00426±0.00066±0.00027 
0.00390±0.00045±0.00030 
0.00282±0.00117±0.00032 
0.00403±0.00079±0.00034 
0.00434±0.00078±0.00036 
0.00450±0.00112±0.00036 
0.00325±0.00080±0.00039 
0.00425±0.00067±0.00040 


0.00064 ±0.00021±0.00006 
0.00057 ±0.00014±0.00007 
0.00026 ±0.00039±0.00008 
0.00058 ±0.00024±0.00008 
0.00067 ±0.00025±0.00009 
0.00080 ±0.00037±0.00009 
0.00044 ±0.00023±0.00009 
0.00066 ±0.00021±0.00010 



Table 9: Mean values and higher moments for the M'^/E'^^^ distributions in the standard, 
15-scheme and p-scheme definitions with statistical and systematic errors. 
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as in 0{ai) 


observable 


1-T 


C 


M^/i^i 


-Dmax 


-Dsum 


average 


fit range 


0.09-0.21 


0.28-0.6 


0.05-0.17 


0.07-0.20 


0.115-0.24 




a,(89.5GeV) 


0.1135 


0.1149 


0.1222 


0.1186 


0.1123 


0.1158 


±A Stat. 
±A sys. exp. 


0.0004 
0.0009 


0.0004 
0.0021 


0.0005 
0.0014 


0.0005 
0.0022 


0.0003 
0.0030 


0.0004 
0.0015 


±A had. 
±A fxji scale 


0.0026 
0.0040 


0.0023 
0.0054 


0.0025 
0.0040 


0.0031 
0.0018 


0.0032 
0.0068 


0.0027 
0.0022 


±Atot 


0.0049 


0.0062 


0.0051 


0.0042 


0.0081 


0.0038 


RMS 


0.0040 


a,(91.2GeV) 


0.1139 


0.1155 


0.1230 


0.1191 


0.1128 


0.1162 


±A Stat. 
±A sys. exp. 


0.0002 
0.0009 


0.0002 
0.0021 


0.0002 
0.0014 


0.0002 
0.0022 


0.0002 
0.0030 


0.0002 
0.0015 


±A had. 
±A /iR scale 


0.0026 
0.0040 


0.0023 
0.0054 


0.0025 
0.0040 


0.0031 
0.0018 


0.0032 
0.0068 


0.0027 
0.0022 


±Atot 


0.0049 


0.0062 


0.0051 


0.0042 


0.0081 


0.0038 


RMS 


0.0042 


a,(93.0GeV) 


0.1128 


0.1140 


0.1208 


0.1179 


0.1117 


0.1152 


±A Stat. 
±A sys. exp. 


0.0004 
0.0009 


0.0004 
0.0021 


0.0005 
0.0014 


0.0005 
0.0022 


0.0003 
0.0030 


0.0003 
0.0015 


±A had. 
±A fiR scale 


0.0026 
0.0040 


0.0023 
0.0054 


0.0025 
0.0040 


0.0031 
0.0018 


0.0032 
0.0068 


0.0027 
0.0022 


±Atot 


0.0049 


0.0062 


0.0051 


0.0042 


0.0081 


0.0038 


RMS 


0.0039 



Table 10: Results of a^ measurements from distributions in 0{a1) 
been reanalyzed for this analysis. 



The data of [10] have 
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a, in 0{ai) 


observable 


1-T 


C 


M^/i?i 


-Dmax 


-Dsum 


average 


fit range 


0.09-0.21 


0.28-0.6 


0.05-0.17 


0.07-0.20 


0.115-0.24 




a,(183GeV) 


0.1075 


0.1053 


0.1053 


0.1065 


0.1024 


0.1059 


±A Stat. 
±A sys. exp. 


0.0042 
0.0022 


0.0039 
0.0023 


0.0047 
0.0024 


0.0044 
0.0026 


0.0034 
0.0012 


0.0037 
0.0020 


±A had. 
±A hr scale 


0.0015 
0.0030 


0.0010 
0.0033 


0.0033 
0.0030 


0.0003 
0.0013 


0.0009 
0.0050 


0.0008 
0.0020 


±Atot 


0.0058 


0.0057 


0.0069 


0.0053 


0.0062 


0.0047 


RMS 


0.0019 


a,(189GeV) 


0.1000 


0.1038 


0.1056 


0.1042 


0.1007 


0.1026 


±A Stat. 
±A sys. exp. 


0.0030 
0.0024 


0.0025 
0.0023 


0.0030 
0.0024 


0.0028 
0.0028 


0.0022 
0.0012 


0.0025 
0.0022 


±A had. 
±A fxji scale 


0.0010 
0.0029 


0.0014 
0.0032 


0.0032 
0.0029 


0.0005 
0.0013 


0.0011 
0.0050 


0.0007 
0.0018 


±Atot 


0.0049 


0.0049 


0.0058 


0.0042 


0.0057 


0.0039 


RMS 


0.0024 


a,(192GeV) 


0.1047 


0.1053 


0.1143 


0.1079 


0.1029 


0.1035 


±A Stat. 
±A sys. exp. 


0.0072 
0.0026 


0.0064 
0.0024 


0.0074 
0.0025 


0.0073 
0.0029 


0.0055 
0.0012 


0.0058 
0.0019 


±A had. 
±A fiR scale 


0.0003 
0.0029 


0.0009 
0.0032 


0.0030 
0.0029 


0.0004 
0.0013 


0.0006 
0.0050 


0.0006 
0.0028 


±Atot 


0.0081 


0.0076 


0.0088 


0.0079 


0.0076 


0.0067 


RMS 


0.0044 


a,(196GeV) 


0.0977 


0.1041 


0.1086 


0.1034 


0.1016 


0.1021 


±A Stat. 
±A sys. exp. 


0.0050 
0.0027 


0.0041 
0.0024 


0.0041 
0.0024 


0.0046 
0.0030 


0.0036 
0.0012 


0.0039 
0.0021 


±A had. 
±A Hfi scale 


0.0011 
0.0029 


0.0009 
0.0032 


0.0038 
0.0029 


0.0003 
0.0013 


0.0005 
0.0049 


0.0005 
0.0023 


±Atot 


0.0065 


0.0058 


0.0072 


0.0057 


0.0062 


0.0050 


RMS 


0.0039 



Table 11: Results of as measurements from distributions in 0{al 
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a, in 0{ai) 


observable 


1-T 


C 


M^/i?i 


-Dmax 


-Dsum 


average 


fit range 


0.09-0.21 


0.28-0.6 


0.05-0.17 


0.07-0.20 


0.115-0.24 




a,(200GeV) 


0.1074 


0.1039 


0.1085 


0.1026 


0.1032 


0.1031 


±A Stat. 
±A sys. exp. 


0.0042 
0.0028 


0.0039 
0.0025 


0.0045 
0.0026 


0.0044 
0.0031 


0.0034 
0.0012 


0.0036 
0.0022 


±A had. 
±A hr scale 


0.0014 
0.0029 


0.0015 
0.0032 


0.0032 
0.0029 


0.0006 
0.0013 


0.0011 
0.0049 


0.0009 
0.0023 


±Atot 


0.0060 


0.0058 


0.0068 


0.0055 


0.0062 


0.0049 


RMS 


0.0027 


a,(202GeV) 


0.1054 


0.1107 


0.1169 


0.1114 


0.1055 


0.1077 


±A Stat. 
±A sys. exp. 


0.0064 
0.0028 


0.0052 
0.0025 


0.0062 
0.0026 


0.0061 
0.0031 


0.0047 
0.0012 


0.0049 
0.0020 


±A had. 
±A fxji scale 


0.0016 
0.0029 


0.0009 
0.0031 


0.0034 
0.0029 


0.0004 
0.0013 


0.0005 
0.0049 


0.0007 
0.0026 


±Atot 


0.0077 


0.0066 


0.0080 


0.0070 


0.0069 


0.0060 


RMS 


0.0048 


a,(205GeV) 


0.0980 


0.0978 


0.1017 


0.1033 


0.0976 


0.1000 


±A Stat. 
±A sys. exp. 


0.0053 
0.0030 


0.0044 
0.0026 


0.0051 
0.0026 


0.0047 
0.0033 


0.0039 
0.0012 


0.0041 
0.0022 


±A had. 
±A fiR scale 


0.0006 
0.0028 


0.0012 
0.0031 


0.0028 
0.0028 


0.0005 
0.0013 


0.0012 
0.0048 


0.0008 
0.0023 


±Atot 


0.0067 


0.0061 


0.0070 


0.0059 


0.0064 


0.0052 


RMS 


0.0026 


a,(207GeV) 


0.1057 


0.1031 


0.1090 


0.1040 


0.1032 


0.1033 


±A Stat. 
±A sys. exp. 


0.0035 
0.0031 


0.0032 
0.0026 


0.0036 
0.0027 


0.0035 
0.0034 


0.0028 
0.0013 


0.0030 
0.0023 


±A had. 
±A Hfi scale 


0.0010 
0.0028 


0.0010 
0.0031 


0.0027 
0.0028 


0.0002 
0.0013 


0.0006 
0.0048 


0.0004 
0.0022 


±Atot 


0.0056 


0.0052 


0.0060 


0.0051 


0.0057 


0.0044 


RMS 


0.0025 



Table 12: Results of a^ measurements from distributions in 0{al 
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a. 


in NLLA 








observable 


1-T 


C 


Mif/^v^is 


-Dmax 


-Dsum 


average 


fit range 


0.09-0.21 


0.16-0.28 


0.05-0.17 


0.07-0.20 


0.115-0.24 




a,(89.5GeV) 


0.1158 


0.1034 


0.1190 


0.1062 


0.1080 


0.1088 


±A Stat. 
±A sys. exp. 


0.0006 
0.0014 


0.0003 
0.0023 


0.0007 
0.0041 


0.0009 
0.0042 


0.0002 
0.0033 


0.0003 
0.0018 


±A had. 
±A X scale 


0.0042 
0.0059 


0.0032 
0.0054 


0.0062 
0.0063 


0.0031 
0.0067 


0.0042 
0.0067 


0.0032 
0.0054 


±Atot 


0.0074 


0.0067 


0.0097 


0.0085 


0.0086 


0.0066 


RMS 


0.0066 


a,(91.2GeV) 


0.1167 


0.1034 


0.1197 


0.10056 


0.1075 


0.1081 


±A Stat. 
±A sys. exp. 


0.0003 
0.0014 


0.0002 
0.0023 


0.0003 
0.0041 


0.0003 
0.0042 


0.0002 
0.0033 


0.0002 
0.0018 


±A had. 
±A X scale 


0.0042 
0.0059 


0.0032 
0.0054 


0.0062 
0.0063 


0.0031 
0.0067 


0.0042 
0.0067 


0.0032 
0.0054 


±Atot 


0.0074 


0.0067 


0.0097 


0.0085 


0.0086 


0.0066 


RMS 


0.0083 


a,(93.0GeV) 


0.1141 


0.1013 


0.1176 


0.1038 


0.1062 


0.1080 


±A Stat. 
d=A sys. exp. 


0.0005 
0.0014 


0.0003 
0.0023 


0.0006 
0.0041 


0.0009 
0.0042 


0.0002 
0.0033 


0.0003 
0.0018 


±A had. 
±A X scale 


0.0042 
0.0059 


0.0032 
0.0054 


0.0062 
0.0063 


0.0031 
0.0067 


0.0042 
0.0067 


0.0032 
0.0054 


±Atot 


0.0074 


0.0067 


0.0097 


0.0085 


0.0086 


0.0066 


RMS 


0.0070 



Table 13: Results of as measurements from distributions in NLLA. The data of [10] have 
been reanalyzed for this analysis. 
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as 


in NLLA 








observable 


1-T 


C 


M^/E-,, 


-Dmax 


-Dsum 


average 


fit range 


0.05-0.09 


0.16-0.28 


0.03-0.05 


0.05-0.07 


0.06-0.115 




a,(183GeV) 


0.0942 


0.0971 


0.1088 


0.1029 


0.1000 


0.1000 


±A Stat. 
±A sys. exp. 


0.0055 
0.0022 


0.0052 
0.0023 


0.0041 
0.0024 


0.0097 
0.0026 


0.0055 
0.0012 


0.0037 
0.0016 


±A had. 
±A X scale 


0.0015 
0.0044 


0.0023 
0.0042 


0.0032 
0.0035 


0.0005 
0.0037 


0.0016 
0.0054 


0.0020 
0.0044 


±Atot 


0.0075 


0.0074 


0.0067 


0.0107 


0.0079 


0.0063 


RMS 


0.0056 


a,(189GeV) 


0.1052 


0.1005 


0.1066 


0.1026 


0.1031 


0.1021 


±A Stat. 
±A sys. exp. 


0.0035 
0.0024 


0.0031 
0.0023 


0.0038 
0.0024 


0.0057 
0.0028 


0.0031 
0.0012 


0.0026 
0.0015 


±A had. 
±A X scale 


0.0025 
0.0044 


0.0025 
0.0041 


0.0031 
0.0035 


0.0005 
0.0037 


0.0019 
0.0053 


0.0021 
0.0041 


±Atot 


0.0066 


0.0062 


0.0065 


0.0074 


0.0065 


0.0055 


RMS 


0.0024 


«,(192GeV) 


0.1058 


0.0906 


0.1007 


0.0890 


0.0998 


0.0998 


±A Stat. 
±A sys. exp. 


0.0101 
0.0026 


0.0090 
0.0024 


0.0110 
0.0025 


0.0099 
0.0029 


0.0082 
0.0012 


0.0071 
0.0016 


±A had. 
±A X scale 


0.0028 
0.0043 


0.0017 
0.0041 


0.0033 
0.0035 


0.0007 
0.0037 


0.0019 
0.0053 


0.0019 
0.0044 


±Atot 


0.0116 


0.0103 


0.0123 


0.0110 


0.0100 


0.0087 


RMS 


0.0071 


a,(196GeV) 


0.1007 


0.0945 


0.0974 


0.1013 


0.0962 


0.0960 


±A Stat. 
±A sys. exp. 


0.0057 
0.0027 


0.0050 
0.0024 


0.0060 
0.0025 


0.0090 
0.0030 


0.0047 
0.0012 


0.0038 
0.0015 


±A had. 
±A X scale 


0.0015 
0.0043 


0.0024 
0.0041 


0.0017 
0.0034 


0.0008 
0.0036 


0.0017 
0.0052 


0.0019 
0.0042 


±Atot 


0.0078 


0.0073 


0.0075 


0.0101 


0.0073 


0.0062 


RMS 


0.0029 



Table 14: Results of a, measurements from distributions in NLLA. 
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as 


in NLLA 








observable 


1-T 


C 


M^/E-,, 


-Dmax 


-Dsum 


average 


fit range 


0.05-0.09 


0.16-0.28 


0.03-0.05 


0.05-0.07 


0.06-0.115 




a,(200GeV) 


0.0981 


0.0864 


0.0974 


0.1012 


0.0953 


0.0914 


±A Stat. 
±A sys. exp. 


0.0058 
0.0028 


0.0047 
0.0025 


0.0057 
0.0026 


0.0085 
0.0031 


0.0045 
0.0014 


0.0038 
0.0017 


±A had. 
±A X scale 


0.0024 
0.0043 


0.0020 
0.0040 


0.0028 
0.0034 


0.0006 
0.0036 


0.0020 
0.0052 


0.0023 
0.0041 


±Atot 


0.0081 


0.0070 


0.0077 


0.0098 


0.0073 


0.0063 


RMS 


0.0056 


a,(202GeV) 


0.1165 


0.1026 


0.1005 


0.1013 


0.1039 


0.1072 


±A Stat. 
±A sys. exp. 


0.0076 
0.0028 


0.0071 
0.0025 


0.0088 
0.0026 


0.0146 
0.0031 


0.0072 
0.0013 


0.0055 
0.0016 


±A had. 
±A X scale 


0.0016 
0.0043 


0.0017 
0.0040 


0.0026 
0.0034 


0.0005 
0.0036 


0.0020 
0.0052 


0.0021 
0.0041 


±Atot 


0.0093 


0.00 87 


0.0101 


0.0154 


0.0092 


0.0073 


RMS 


0.0066 


«,(205GeV) 


0.0928 


0.0970 


0.1032 


0.1039 


0.1036 


0.0996 


±A Stat. 
±A sys. exp. 


0.0056 
0.0030 


0.0050 
0.0026 


0.0062 
0.0026 


0.0098 
0.0031 


0.0049 
0.0013 


0.0039 
0.0015 


±A had. 
±A X scale 


0.0029 
0.0042 


0.0023 
0.0040 


0.0024 
0.0034 


0.0005 
0.0036 


0.0017 
0.0051 


0.0022 
0.0041 


±Atot 


0.0081 


0.0073 


0.0079 


0.0109 


0.0074 


0.0063 


RMS 


0.0050 


a,(207GeV) 


0.1054 


0.0935 


0.1012 


0.0972 


0.0975 


0.0976 


±A Stat. 
±A sys. exp. 


0.0043 
0.0031 


0.0038 
0.0026 


0.0047 
0.0027 


0.0066 
0.0034 


0.0036 
0.0013 


0.0030 
0.0016 


±A had. 
±A X scale 


0.0019 
0.0042 


0.0018 
0.0040 


0.0025 
0.0034 


0.0005 
0.0036 


0.0016 
0.0051 


0.0019 
0.0040 


±Atot 


0.0070 


0.0064 


0.0069 


0.0083 


0.0066 


0.0056 


RMS 


0.0045 



Table 15: Results of a, measurements from distributions in NLLA. 



29 



as in C>(a,")+NLLA (logR) 


observable 


1-T 


C 


M^/i54 


-Dmax 


-Dsum 


average 


fit range 


0.09-0.21 


0.16-0.6 


0.05-0.17 


0.07-0.20 


0.115-0.24 




a,(89.5GeV) 


0.1257 


0.1211 


0.1226 


0.1155 


0.1249 


0.1220 


±A Stat. 
±A sys. exp. 


0.0002 
0.0021 


0.0002 
0.0022 


0.0004 
0.0023 


0.0004 
0.0021 


0.0001 
0.0019 


0.0004 
0.0020 


±A had. 
±A X scale 


0.0031 
0.0060 


0.0021 
0.0057 


0.0019 
0.0048 


0.0023 
0.0051 


0.0029 
0.0073 


0.0020 
0.0048 


±Atot 


0.0071 


0.0065 


0.0057 


0.0060 


0.0081 


0.0056 


RMS 


0.0040 


a,(91.2GeV) 


0.1256 


0.1211 


0.1230 


0.1156 


0.1250 


0.1219 


±A Stat. 
±A sys. exp. 


0.0002 
0.0021 


0.0002 
0.0022 


0.0002 
0.0023 


0.0001 
0.0021 


0.0001 
0.0019 


0.0002 
0.0020 


±A had. 
±A X scale 


0.0031 
0.0060 


0.0021 
0.0057 


0.0019 
0.0048 


0.0023 
0.0051 


0.0029 
0.0073 


0.0020 
0.0048 


±Atot 


0.0071 


0.0065 


0.0057 


0.0060 


0.0081 


0.0056 


RMS 


0.0040 


a,(93.0GeV) 


0.1257 


0.1211 


0.1208 


0.1144 


0.1235 


0.1222 


±A Stat. 
±A sys. exp. 


0.0002 
0.0021 


0.0002 
0.0022 


0.0004 
0.0023 


0.0004 
0.0021 


0.0001 
0.0019 


0.0004 
0.0020 


±A had. 
±A X scale 


0.0031 
0.0060 


0.0021 
0.0057 


0.0019 
0.0048 


0.0023 
0.0051 


0.0029 
0.0073 


0.0020 
0.0048 


±Atot 


0.0071 


0.0065 


0.0057 


0.0060 


0.0081 


0.0056 


RMS 


0.0042 



Table 16: Results of a^ measurements from distributions in (9(a^)+NLLA. The data 
of [10] have been reanalyzed for this analysis. 
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as in C(a^)+NLLA (logR) 


observable 


1-T 


C 


M^E'L 


-Dmax 


-Dsum 


average 


fit range 


0.05-0.21 


0.16-0.6 


0.03-0.17 


0.05-0.20 


0.06-0.24 




a,(183GeV) 


0.1072 


0.1111 


0.1094 


0.1049 


0.1119 


0.1081 


±A Stat. 
±A sys. exp. 


0.0043 
0.0022 


0.0036 
0.0023 


0.0037 
0.0024 


0.0037 
0.0026 


0.0039 
0.0012 


0.0034 
0.0021 


±A had. 
±A X scale 


0.0017 
0.0044 


0.0017 
0.0042 


0.0031 
0.0035 


0.0005 
0.0037 


0.0016 
0.0054 


0.0011 
0.0037 


±Atot 


0.0068 


0.0062 


0.0064 


0.0059 


0.0070 


0.0056 


RMS 


0.0029 


a,(189GeV) 


0.1105 


0.1101 


0.1079 


0.1032 


0.1122 


0.1067 


±A Stat. 
±A sys. exp. 


0.0027 
0.0024 


0.0023 
0.0023 


0.0027 
0.0024 


0.0024 
0.0028 


0.0021 
0.0012 


0.0022 
0.0020 


±A had. 
±A X scale 


0.0020 
0.0044 


0.0020 
0.0041 


0.0031 
0.0035 


0.0005 
0.0037 


0.0016 
0.0053 


0.0010 
0.0038 


±Atot 


0.0060 


0.0056 


0.0059 


0.0052 


0.0060 


0.0049 


RMS 


0.0035 


a,(192GeV) 


0.1124 


0.1083 


0.1082 


0.1051 


0.1139 


0.1096 


±A Stat. 
±A sys. exp. 


0.0068 
0.0026 


0.0060 
0.0024 


0.0063 
0.0025 


0.0060 
0.0029 


0.0052 
0.0012 


0.0051 
0.0018 


±A had. 
±A X scale 


0.0017 
0.0043 


0.0014 
0.0041 


0.0032 
0.0035 


0.0007 
0.0037 


0.0019 
0.0053 


0.0015 
0.0040 


±Atot 


0.0086 


0.0078 


0.0083 


0.0077 


0.0078 


0.0069 


RMS 


0.0035 


«,(196GeV) 


0.1045 


0.1079 


0.1060 


0.1024 


0.1107 


0.1068 


±A Stat. 
±A sys. exp. 


0.0042 
0.0027 


0.0038 
0.0024 


0.0040 
0.0025 


0.0038 
0.0030 


0.0032 
0.0012 


0.0033 
0.0019 


±A had. 
±A X scale 


0.0013 
0.0043 


0.0013 
0.0041 


0.0029 
0.0034 


0.0008 
0.0036 


0.0017 
0.0052 


0.0014 
0.0038 


±Atot 


0.0067 


0.0062 


0.0065 


0.0061 


0.0065 


0.0055 


RMS 


0.0036 



Table 17: Results of ag measurements from distributions in 0{af)+NLLA. 
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as in C(a^)+NLLA (logR) 


observable 


1-T 


C 


M^E'L 


-Dmax 


-Dsum 


average 


fit range 


0.05-0.21 


0.16-0.6 


0.03-0.17 


0.05-0.20 


0.06-0.24 




a,(200GcV) 


0.1095 


0.1044 


0.1045 


0.1021 


0.1101 


0.1044 


±A Stat. 
±A sys. exp. 


0.0041 
0.0028 


0.0036 
0.0025 


0.0034 
0.0026 


0.0036 
0.0031 


0.0032 
0.0012 


0.0031 
0.0020 


±A had. 
±A X scale 


0.0023 
0.0043 


0.0022 
0.0040 


0.0030 
0.0034 


0.0006 
0.0036 


0.0020 
0.0052 


0.0015 
0.0037 


±Atot 


0.0071 


0.0063 


0.0062 


0.0060 


0.0066 


0.0054 


RMS 


0.0036 


a,(202GeV) 


0.1188 


0.1160 


0.1105 


0.1103 


0.1194 


0.1141 


±A Stat. 
±A sys. exp. 


0.0058 
0.0028 


0.0050 
0.0025 


0.0055 
0.0026 


0.0053 
0.0031 


0.0047 
0.0012 


0.0046 
0.0019 


±A had. 
±A X scale 


0.0018 
0.0043 


0.0015 
0.0040 


0.0031 
0.0034 


0.0005 
0.0036 


0.0016 
0.0052 


0.0014 
0.0038 


±Atot 


0.0081 


0.0070 


0.0074 


0.0071 


0.0074 


0.0064 


RMS 


0.0044 


a,(205GeV) 


0.1023 


0.1064 


0.1041 


0.1031 


0.1109 


0.1071 


±A Stat. 
±A sys. exp. 


0.0045 
0.0030 


0.0037 
0.0026 


0.0041 
0.0026 


0.0038 
0.0033 


0.0033 
0.0012 


0.0033 
0.0019 


±A had. 
±A X scale 


0.0021 
0.0042 


0.0017 
0.0040 


0.0027 
0.0034 


0.0005 
0.0036 


0.0017 
0.0051 


0.0012 
0.0037 


±Atot 


0.0072 


0.0063 


0.0065 


0.0061 


0.0064 


0.0055 


RMS 


0.0035 


«,(207GeV) 


0.1118 


0.1074 


0.1058 


0.1021 


0.1134 


0.1061 


±A Stat. 
±A sys. exp. 


0.0036 
0.0031 


0.0026 
0.0026 


0.0031 
0.0027 


0.0029 
0.0034 


0.0025 
0.0013 


0.0026 
0.0020 


±A had. 
±A X scale 


0.0016 
0.0042 


0.0014 
0.0040 


0.0027 
0.0034 


0.0005 
0.0036 


0.0016 
0.0051 


0.0012 
0.0037 


±Atot 


0.0064 


0.0056 


0.0060 


0.0058 


0.0060 


0.0051 


RMS 


0.0050 



Table 18: Results of ag measurements from distributions in 0{af)+NLLA. 
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Us from mean values with power corrections 


Observable 


(1-T) 


(C) 


{M-^/E-D 


(-Bmax) 


(Bsum) 


average 


a,(45GeV) 


0.1341 


0.1418 


0.1268 


0.1358 


0.1421 


0.1370 


±A Stat. 
±A sys. exp. 


0.0154 
0.0036 


0.0165 
0.0037 


0.0140 
0.0026 


0.0127 
0.0005 


0.0091 
0.0005 


0.0043 
0.0023 


±A /i/j scale 
±A /i/ scale 


0.0072 
0.0052 


0.0055 
0.0038 


0.0041 
0.0035 


0.0062 
0.0005 


0.0052 
0.0001 


0.0051 
0.0016 


±Atot 


0.0181 


0.0181 


0.0152 


0.0142 


0.0105 


0.0073 


RMS 


0.0063 


a,(66GeV) 


0.1159 


0.1252 


0.1140 


0.1189 


0.1265 


0.1251 


±A Stat. 
±A sys. exp. 


0.0080 
0.0023 


0.0081 
0.0028 


0.0076 
0.0032 


0.0070 
0.0009 


0.0055 
0.0006 


0.0043 
0.0010 


±A ^R scale 
±A Hi scale 


0.0059 
0.0054 


0.0045 
0.0049 


0.0034 
0.0030 


0.0051 
0.0012 


0.0043 
0.0009 


0.0041 
0.0019 


±Atot 


0.0116 


0.0109 


0.0094 


0.0088 


0.0070 


0.0064 


RMS 


0.0055 


03(76 GeV) 


0.1302 


0.1388 


0.1235 


0.1244 


0.1313 


0.1255 


±A Stat. 
±A sys. exp. 


0.0074 
0.0032 


0.0077 
0.0040 


0.0072 
0.0034 


0.0070 
0.0020 


0.0053 
0.0016 


0.0039 
0.0023 


±A Hfi scale 
±A III scale 


0.0055 
0.0024 


0.0042 
0.0013 


0.0032 
0.0018 


0.0048 
0.0006 


0.0040 
0.0004 


0.0038 
0.0008 


±Atot 


0.0101 


0.0098 


0.0087 


0.0088 


0.0068 


0.0060 


RMS 


0.0062 


a,(89GeV) 


0.1189 


0.1263 


0.1152 


0.1163 


0.1235 


0.1177 


±A Stat. 
±A sys. exp. 


0.0004 
0.0015 


0.0004 
0.0013 


0.0004 
0.0019 


0.0003 
0.0003 


0.0002 
0.0004 


0.0003 
0.0009 


±A Hn scale 
±A yU/ scale 


0.0052 
0.0031 


0.0039 
0.0027 


0.0029 
0.0019 


0.0044 
0.0009 


0.0037 
0.0007 


0.0031 
0.0010 


±Atot 


0.0063 


0.0050 


0.0040 


0.0046 


0.0038 


0.0034 


RMS 


0.0048 


as(91.2GeV) 


0.1193 


0.1270 


0.1153 


0.1167 


0.1238 


0.1176 


±A Stat. 
±A sys. exp. 


0.0002 
0.0015 


0.0001 
0.0013 


0.0001 
0.0019 


0.0001 
0.0008 


0.0001 
0.0007 


0.0001 
0.0011 


±A hr scale 
±A Hi scale 


0.0051 
0.0030 


0.0039 
0.0025 


0.0029 
0.0019 


0.0044 
0.0008 


0.0037 
0.0007 


0.0031 
0.0010 


±Atot 


0.0061 


0.0048 


0.0039 


0.0046 


0.0038 


0.0034 


RMS 


0.0049 



Table 19: Results of ag measurements from mean values with power corrections. The 
data of [10] and [33] have been reanalyzed for this analysis 
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Us from mean values with power corrections 


Observable 


(1-T) 


(C) 


m/m,:) 


(-Bmax) 


X-Dsum/ 


average 


^.(gSGeV) 


0.1182 


0.1255 


0.1145 


0.1157 


0.1226 


0.1171 


±A Stat. 
±A sys. exp. 


0.0004 
0.0013 


0.0003 
0.0011 


0.0004 
0.0017 


0.0003 
0.0006 


0.0002 
0.0004 


0.0003 
0.0009 


±A ^R scale 
±A yU/ scale 


0.0051 
0.0030 


0.0039 
0.0026 


0.0029 
0.0019 


0.0044 
0.0009 


0.0037 
0.0007 


0.0031 
0.0010 


±Atot 


0.0061 


0.0048 


0.0039 


0.0045 


0.0038 


0.0033 


RMS 


0.0047 


as(133GeV) 


0.1158 


0.1203 


0.1120 


0.1109 


0.1163 


0.1150 


±A Stat. 
±A sys. exp. 


0.0047 
0.0010 


0.0039 
0.0011 


0.0048 
0.0009 


0.0038 
0.0016 


0.0026 
0.0011 


0.0023 
0.0011 


±A yU/j scale 
±A ^i scale 


0.0043 
0.0020 


0.0033 
0.0019 


0.0025 
0.0013 


0.0037 
0.0007 


0.0031 
0.0006 


0.0030 
0.0010 


±Atot 


0.0067 


0.0056 


0.0056 


0.0056 


0.0043 


0.0041 


RMS 


0.0037 


a,(161GeV) 


0.1037 


0.1131 


0.1020 


0.1053 


0.1083 


0.1068 


±A Stat. 
±A sys. exp. 


0.0069 
0.0052 


0.0067 
0.0066 


0.0068 
0.0060 


0.0053 
0.0038 


0.0038 
0.0023 


0.0031 
0.0026 


±A iiR scale 
±A yU/ scale 


0.0040 
0.0023 


0.0031 
0.0020 


0.0023 
0.0014 


0.0034 
0.0007 


0.0029 
0.0008 


0.0028 
0.0010 


±Atot 


0.0098 


0.0101 


0.0094 


0.0074 


0.0053 


0.0050 


RMS 


0.0044 


a,(172GeV) 


0.1140 


0.1164 


0.1107 


0.1119 


0.1167 


0.1130 


±A Stat. 
±A sys. exp. 


0.0112 
0.0021 


0.0090 
0.0036 


0.0106 
0.0011 


0.0080 
0.0020 


0.0069 
0.0048 


0.0061 
0.0031 


±A yU,R scale 
±A yU/ scale 


0.0039 
0.0015 


0.0030 
0.0016 


0.0022 
0.0010 


0.0034 
0.0004 


0.0028 
0.0004 


0.0031 
0.0011 


±Atot 


0.0121 


0.0103 


0.0109 


0.0089 


0.0089 


0.0076 


RMS 


0.0027 



Table 20: Results of a^ measurements from mean values with power corrections, 
data of [10] have been reanalyzed for this analysis 



The 
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Us from mean values with power corrections 


observable 


(1-T) 


(C) 


{M-^/Ei:) 


(-Bmax) 


(Bsum) 


average 


as(183GeV) 


0.1156 


0.1172 


0.1067 


0.1059 


0.1121 


0.1121 


±A Stat. 
±A sys. exp. 


0.0047 
0.0045 


0.0038 
0.0061 


0.0058 
0.0032 


0.0042 
0.0044 


0.0026 
0.0011 


0.0022 
0.0013 


±A hr scale 
±A III scale 


0.0038 
0.0013 


0.0029 
0.0014 


0.0022 
0.0010 


0.0033 
0.0005 


0.0027 
0.0005 


0.0029 
0.0005 


±Atot 


0.0077 


0.0079 


0.0071 


0.0069 


0.0040 


0.0039 


RMS 


0.0051 


as(189GeV) 


0.1092 


0.1185 


0.1022 


0.1039 


0.1085 


0.1082 


±A Stat. 
±A sys. exp. 


0.0032 
0.0044 


0.0024 
0.0061 


0.0040 
0.0032 


0.0028 
0.0044 


0.0017 
0.0011 


0.0017 
0.0012 


±A yU/j scale 
±A yU/ scale 


0.0037 
0.0016 


0.0029 
0.0013 


0.0021 
0.0011 


0.0032 
0.0006 


0.0027 
0.0006 


0.0027 
0.0006 


±Atot 


0.0068 


0.0073 


0.0057 


0.0061 


0.0035 


0.0034 


RMS 


0.0064 


a,(192GeV) 


0.0984 


0.1056 


0.0947 


0.0980 


0.1055 


0.1092 


±A Stat. 
±A sys. exp. 


0.0093 
0.0040 


0.0066 
0.0057 


0.0127 
0.0030 


0.0074 
0.0042 


0.0045 
0.0010 


0.0028 
0.0015 


±A yU/j scale 
±A yU/ scale 


0.0037 
0.0021 


0.0028 
0.0021 


0.0021 
0.0013 


0.0032 
0.0007 


0.0027 
0.0007 


0.0030 
0.0009 


±Atot 


0.0110 


0.0094 


0.0133 


0.0091 


0.0054 


0.0044 


RMS 


0.0049 


a,(196GeV) 


0.1165 


0.1124 


0.1036 


0.1072 


0.1100 


0.1092 


±A Stat. 
±A sys. exp. 


0.0057 
0.0045 


0.0044 
0.0059 


0.0070 
0.0032 


0.0048 
0.0044 


0.0030 
0.0011 


0.0023 
0.0014 


±A hr scale 
±A yU/ scale 


0.0037 
0.0012 


0.0028 
0.0016 


0.0021 
0.0010 


0.0032 
0.0005 


0.0027 
0.0005 


0.0029 
0.0006 


±Atot 


0.0082 


0.0080 


0.0081 


0.0072 


0.0042 


0.0039 


RMS 


0.0049 



Table 21: Results of a^ measurements from mean values with power corrections. 
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as from mean values with power corrections 


observable 


(1-T) 


(C) 


{Mi/Ei:) 


(-Bmax) 


(Bsum) 


average 


as(200GeV) 


0.1068 


0.1104 


0.1062 


0.1056 


0.1096 


0.1105 


±A Stat. 
±A sys. exp. 


0.0057 
0.0042 


0.0043 
0.0058 


0.0069 
0.0032 


0.0048 
0.0044 


0.0030 
0.0011 


0.0023 
0.0014 


±A hr scale 
±A III scale 


0.0037 
0.0016 


0.0028 
0.0017 


0.0021 
0.0009 


0.0032 
0.0005 


0.0026 
0.0005 


0.0028 
0.0006 


±Atot 


0.0082 


0.0079 


0.0079 


0.0072 


0.0042 


0.0039 


RMS 


0.0022 


as(202GeV) 


0.0946 


0.1066 


0.1085 


0.1070 


0.1127 


0.1185 


±A Stat. 
±A sys. exp. 


0.0097 
0.0040 


0.0063 
0.0057 


0.0101 
0.0032 


0.0067 
0.0044 


0.0040 
0.0011 


0.0023 
0.0015 


±A yU/j scale 
±A yU/ scale 


0.0036 
0.0022 


0.0028 
0.0019 


0.0021 
0.0009 


0.0031 
0.0005 


0.0026 
0.0004 


0.0029 
0.0008 


±Atot 


0.0114 


0.0091 


0.0108 


0.0086 


0.0049 


0.0041 


RMS 


0.0068 


as(205GeV) 


0.0877 


0.0879 


0.0942 


0.0991 


0.1003 


0.1042 


±A Stat. 
±A sys. exp. 


0.0072 
0.0037 


0.0054 
0.0052 


0.0078 
0.0030 


0.0050 
0.0042 


0.0033 
0.0011 


0.0023 
0.0013 


±A yU/j scale 
±A yU/ scale 


0.0036 
0.0025 


0.0028 
0.0030 


0.0021 
0.0012 


0.0031 
0.0006 


0.0026 
0.0008 


0.0028 
0.0010 


±Atot 


0.0092 


0.0085 


0.0087 


0.0072 


0.0044 


0.0040 


RMS 


0.0060 


a,(207GeV) 


0.1063 


0.1077 


0.1049 


0.1055 


0.1070 


0.1072 


±A Stat. 
±A sys. exp. 


0.0047 
0.0042 


0.0036 
0.0057 


0.0055 
0.0032 


0.0040 
0.0044 


0.0024 
0.0010 


0.0020 
0.0012 


±A hr scale 
±A yU/ scale 


0.0036 
0.0015 


0.0028 
0.0018 


0.0020 
0.0009 


0.0031 
0.0005 


0.0026 
0.0006 


0.0027 
0.0006 


±Atot 


0.0074 


0.0075 


0.0068 


0.0067 


0.0037 


0.0036 


RMS 


0.0011 



Table 22: Results of a^ measurements from mean values with power corrections. 
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a( 192 GeV) from distributions in 0(a 2) 
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mX' 
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B 
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66GeV (0.01) 

76GeV (-0.01) 

89GeV (0.25) 

91 GeV ( 0.17) 

93GeV (0.28) 

133 GeV ( 0.02) 

161 GeV (-0.02) 

172 GeV (-0.01) 

183 GeV ( 0.06) 

189 GeV ( 0.19) 

192 GeV (-0.01) 

196 GeV ( 0.03) 

200 GeV ( 0.03) 

202 GeV ( 0.01) 

205 GeV (-0.05) ^ 
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DELPHI 
205 GeV 

• data 
Jetset 7.4 PS 
Herwig 5.8 
Ariadne 4.08 




B 



sum 



Forward Chamber A 



Forward RICH 



Barrel Muon Chambers 



Forward Chamber B 



Forward EM Calorimeter 

Forward Hadron Calorimeter 

Forward Hodoscope 

Forward Muon Chambers 

Surround Muon Chambers 



DELPHI 



Barrel Hadron Calorimeter 



Scintillators 




Superconducting Coil 

h Density Projection Chamber 
Outer Detector 

Barrel RICH 

Small Angle Tile Calorimeter 



Quadrupole 



Very Small Angle Tagger 



Beam Pipe 
Vertex Detector 



Inner Detector 



Time Projection Chamber 




M, . 2/E . 2 
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